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Abstract 

We investigate a U(l) chiral gauge model in 4+1 dimensions formulated 

on the lattice via the domain-wall method. We calculate an effective action 

for smooth background gauge fields at a fermion one loop level. From this 

calculation we discuss properties of the resulting 4 dimensional theory, such 

as gauge invariance of 2 point functions, gauge anomalies and an anomaly in 

the fermion number current. 
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I. INTRODUCTION 



The standard model is very successful to explain many aspects of electro-weak interac- 
tions. However these successes come mainly from perturbative analysis, and physics at the 
breaking scale, for example, the mass of the Higgs particle and the baryon number violation 
etc, can not be predicted. In order to predict them, we need to study the standard model 
non-perturbatively - especially using the technique of lattice gauge theories. 

A main problem for studying the standard model on the lattice is the difficulty to define 
lattice chiral gauge theories due to the fermion doubling phenomenon which can be 

easily seen in the fermion propagator on the lattice; 



This propagator has poles at p M = (tt, 0, 0, 0), etc as well as at p M = (0, 0, 0, 0). Therefore a 
naively discretized lattice fermion field yields 2 d fermion modes, half of one chirality and half 
of the other, so that the theory is no more chiral and therefore can not be used to construct 
the standard model on the lattice. Several lattice approaches have been proposed to define 
chiral gauge theories, but so far none of them have been proven to work successfully. 

Recently Kaplan has proposed a new approach || to this problem. He suggested that 
it may be possible to simulate the behavior of massless chiral fermions in 2k dimensions by 
a lattice theory of massive fermions in 2k+l dimensions if the fermion mass has a shape 
of a domain wall in the 2k+l-th dimension. He showed for the weak gauge coupling limit 
that the massless chiral states arise as zero-modes bound to the 2k-dimensional domain wall 
while all doublers can be given large gauge invariant masses. If the chiral fermion content 
that appears on the domain wall is anomalous the 2k-dimensional gauge current flows off 
the wall into the extra dimension so that the theory can not be 2k-dimensional. Therefore 
he argued that this approach possibly simulates the 2k-dimensional chiral fermions only for 
anomaly-free cases. 

His idea, called a domain-wall fermion method, was tested for smooth external gauge 
fields. It has been shown both numerically and analytically [f| that in the case of the 
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chiral Schwinger model the anomaly in the gauge current is cancelled on the wall among 
three fermions of charge 3, 4, and 5. The Chern-Simons current was also evaluated in Ref. 
|6|||: It is shown that the 2k+l-th component of the current is non-zero in the positive mass 
region and zero in the negative mass region, so that the derivative of the current cancels the 
2k-dimensional gauge anomaly on the wall, as was argued in Ref. |3|,0] []. 

Results above provide positive indications that the domain-wall fermion method may 
work as a lattice regularization for chiral gauge theories. There exists, however, two remain- 
ing problems to be considered in this approach. One of the problems is the fate of the chiral 
zero mode on the domain wall: Since the original 2k+l-dimensional model is vector-like, 
there always exists an anti-chiral mode, localized on an anti-domain wall formed by period- 
icity of the extra dimension. If the chiral mode and the anti-chiral mode are paired into a 
Dirac mode, this approach fails to simulate chiral gauge theories. Without dynamical gauge 
fields, the overlap between the chiral mode and the anti-chiral mode is suppressed as 0(e~ L ) 
where L is the size of the extra dimension. If gauge fields become dynamical, the overlap 
depends on the gauge coupling. It was found |TT| , |T2"[| that the chiral mode disappears and 



the model becomes vector-like in the strong gauge coupling limit of the extra dimension. 
Recently this problem has been investigated at the intermediate coupling region via the 
numerical simulation for a 2+1 dimensional U(l) model, but no definite conclusion on the 



existence of the chiral zero mode can be obtained in the symmetric phase |L3 . 

The other problem is related to a structure of an effective action for smooth back-ground 
gauge fields at the fermion 1-loop level: The perturbative evaluation for the 2+1 dimen- 
sional model found that, if gauge fields depend on coordinates of the extra dimension, 
the effective action contains the longitudinal component as well as parity-odd terms, and 
that this longitudinal component, which breaks gauge invariance, remains nonzero even for 



1 Recently the anomaly is also calculated in the continuum version of the domain-wall fermion 
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anomaly-free cases. The gauge non-invariant parity-even term seems absent [14|,l5|l in two 



modifications of the Kaplan's original domain wall fermion, wave-guide model ||T6| , |1~7|| and 
overlap formula p!8|JT9|| , whose gauge fields do not depend on coordinates of the extra di- 
mension. These two modifications, however, suffers from the first problem: the chiral zero 
mode on the domain wall seems to disappear in the presence of the dynamical gauge fields 



In this paper, following the previous calculation in 2+1 dimensions ||, we have carried 
out a detailed perturbative calculation of the original domain wall fermion formulation in 
4+1 dimensions for smooth background gauge fields, in order to investigate the structure 
of the effective action in higher dimensions. In sect, we briefly summarize the lattice 
perturbation theory for the domain wall method with the periodic boundary condition ||. 
In sect. |T|, we evaluate the 2-point function and anomaly in 4+1 dimensions at a fermion 
one-loop order. We find that the effective action for the 4+1 dimensional theory has the 
similar structure to that for the 2+1 dimensional theory: there appear not only parity-odd 
terms such as the gauge anomaly and the Chern-Siomns term but also parity even terms such 
as the mass term and the Lorentz non-covariant term. Therefore the gauge non-invariant 
terms remin non-zero for anomaly free cases also in 4 dimensions. In sect. [TV], we comment 
on an anomaly of the fermion number current in 4+1 dimensions. Finally we give our 
conclusions in sect. |V|. 

II. FORMULATION 

In this section, we briefly summarize a formulation of the domain wall method and a 
set-up of lattice perturbation theories. In particular, we explicitly give a fermion propagator, 
vertex functions, and the Ward-Takahashi identity. 
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A. Lattice Action 



We consider a vector gauge theory in D=2k+1 dimensions with a domain wall mass term. 



For later convenience we use the notation of Ref. 18 , where the fermionic action is written 



in terms of a 2k dimensional theory with infinitely many flavors. Our action is denoted as 

S = S G + S F . (2) 
The action for gauge field Sq is given by 

S g = P E E Re {Tr[£V(n, S )]} 

n,fi>v s 

+ ^EE Re {Tr[^D(n,s)]} (3) 

71, fl S 

where ft, v run from 1 to 2k, n is a point on a 2k-dimensional lattice and s a coordinate 
in the extra dimension, (3 is the inverse gauge coupling for plaquettes U^ v and (3d that for 
plaquettes U^d- The fermionic part of the action Sf is given by 

Sf = \ E E ^s{n)-f^[U s ^{n)tp s (n + fi) 

n,fi s 

+ ^2J2Mn)[M P R + MlP L } st i, t {n) (4) 

n s,t 

+ ^J2J2^s(n)[U s ^(n)i) s {n + fi) 

n,fi s 

+ u l,M - v)*Ps(n -fi)- 2ip s (n)} 
where s, t are considered as flavor indices, Pr/l — (1 ± 72fc+i)/2, 

(M )st = f/ 8 ,D(n)(5 s+ i,t - a(s)5 si 

(M % = E7J_i,dH^-i,< " (5) 

and U Sttl (n), U s ^{n) are link variables for gauge fields. We consider the above model with a 
periodic boundary in the extra dimension, so that s, t run from — L to L — 1, and we take 
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a(s) = 1 - m [sign(s + -) • sign(L - s - -)] 

1 — m , — \ < s < L — ~ 
1 + m , — L — | < s < — | 



(6) 



for — L < s < L. It is easy to see |18[ that Sf above is identical to the domain-wall fermion 
action in D=2k+1 dimensions with the Wilson parameter r = 1. In fact the second term 
in eq. (HI) can be rewritten as 



(7) 



with M(s) = m [sign(s + 1/2) ■ sign(L — s — 1/2)] . It is easy to see that one chiral zero mode 
appears on the domain wall around s — if and only if < m < 2. 



B. Fermion propagator and Feynman rules 

The fermion propagator in 2k-dimensional momentum space and in real D-th space has 
been obtained in Ref. [H|5[] for large L: 



where 



G L (p) 



St 



S F (p) st = -[[(^W + M)G L (p)] st P L 
+[(*E7^ + M*)G R (p)UP R ] 

[s,t> 0) 

Be -a+\s-t\ + ( Al _ B y-a+{s+t) + ( Ar _ B y~a+{2L-s~t)^ 

[s>0, t < 0) 

[s < 0, t > 0) 

A L e a - s ~ a+t + y 4 R e~ a -( L+s )~ a+(7 '~ t ), 
[s,t< 0) 

Ce -a-\»-t\ + (A L - C)e a - (s+i) + (A R - C)e- a - {2L+s+t \ 



(8) 



(9) 



G R (p) 



St 



with 



(s,t > -1) 

Be -a+\s-t\ + _ B y- a+ {s+t+2) 

+ {A L - B )e- a +^ 2L - s - t - 2 \ 

(s >-l,t< -1) 

A Re -a+(s+l)+a-(t+l) _|_ ^ ie -a+(L-s-l)-a_(L+t+l) 

(s < -l, * > -1) 

^ e a-(s+l)-a + (t+l) + ^4 ie -Q_(L+ S +l)-a + (L-t-l)^ 

(s,f < -1) 

C e -a_| a -t| + ( Ar _ C ) e a_( S+t +2) 

+ (A L -C)e^-( 2L+s+ '+ 2 ), 



, V(p) 
a± = 1 — =F m 



1 l+p 2 

a± = arccosh[-(a.|- H )] > 0, 

2 Q-j- 



B = 
C = 



1 



1 

a_e a - — a_|_e~ a + 
1 



2a + sinho; + ' 
1 

2a_sinho;_ 



(10) 



(11) 
(12) 

(13) 

(14) 

(15) 

(16) 



From the form of A L , A R , B and C, it is easy to see that singularities occur only in A L 
at p = 0; 

m (4 - ml) 



A, 



Ap 2 a? 



p^O. 



(17) 



Therefore the Gl part of the propagator describes one massless right-handed fermion around 
s,t — 0, which corresponds to the zero mode on the domain wall. It is also noted that the 
Gr part describes one massless left-handed fermion around \s\, \t\ = L, which corresponds 
to the anti-zero mode, due to the periodic boundary condition in the extra dimension. 



Now we write down the lattice Feynman rules relevant for a fermion one-loop calculation, 
which will be performed in the next section. We first choose the axial gauge fixing^ : 
U s d = 1. Although the full gauge symmetries in D dimensions are lost, the theory is 



still invariant under s-independent gauge transformations |zl|| . Therefore the gauge current 
Jfi( x ) — * s conserved. We consider the limit of small 2k-dimensional gauge 

coupling, and take 

U s>II (n) = exp[iagA^(s, n + /x/2)] (18) 

where a is the lattice spacing, and g oc is the gauge coupling constant whose mass 

dimension is 2 — D/2 (mass dimension of the gauge fields is D/2 — 1). We consider 
Feynman rules in momentum space for the physical 2k dimensions but in real space for the 
extra dimension. There are three relevant points for later calculations. 

The fermion propagator SV(p)st has been given in eq.(j§) with eqs.(|| p!0|) . 

The fermion vertex coupled to a single gauge field is given by 

agd fl [S F 1 (k)] st = i cos (ak^)j ^5 Stt + sin(a^)5 Sit (19) 

q + p 

where k = and the fermion vertex with two gauge fields is 

2 

- a^dl[S F \ q -±P)] ss . (20) 

From the periodicity of eq. (|T9"D , the fermion vertex with 2n + 1 gauge fields is proportional to 
^[S , i ? 1 ( 2 ^ £ )] st and the fermion vertex with 2n gauge fields is proportional to ^[S , ^ 1 ( 2 ^ £ )] st . 



2 To choose the axial gauge fixing in the periodic boundary condition, gauge field configurations 
should satisfy a constraint that the Polyakov loop in the extra dimension is equal to unity. To 
achieve this constraint, we should put a delta function of the constraint, so that the other gauge 
coupling g s oc is not necessarily small and can be made arbitrary large, or we should also 

take the weak coupling limit of g s . 
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The fermion propagator satisfies the Ward-Takahashi identity on the lattice. The identity 
is given by 

f^2sin(j ) ,/2)d,[S F 1 (q + P /2)] s , t = [S F \q + p)] s , t 

" [S?\q)Ut. (21) 



III. U(l) CHIRAL GAUGE THEORY IN 4+1 DIMENSIONS 

In this section we investigate an U(l) chiral gauge theory in 4+1 dimensions. In 4 dimen- 
sions, it is shown with the power counting that the n-point functions which has divergent 
diagrams are 2-, 3- and 4-point functions. In the following two subsections we calculate in 
detail the 2-point function of the gauge field Ii^ v (j)) and the parity odd part of the 3-point 
function ( the gauge anomaly ). 

A. Calculation of 2-point function 

First of all, the effective action with two external gauge fields is denoted by 

s ?ff = ~C E M^p)Mt, -vW u {v)st 

P,3,t 

= ~E A,(s,p)A u (t, -p)[ITt 2 ) + Ilf (22) 

P,3,t 



where 



[^Tst = f_'* a (^^{[^(ff + f ) ■ S F {q + V )Ud v S- F \q + §) ■ S F {q)]ts) X a\ (23) 

[iffe = -W^ f /a f ^AdlS F \q) ■ S F (q)] ss x a 2 . (24) 

This integral has the similar form as in the 2+1 dimensions case, but has the divergence 
of order a~ 2 . To separate the would-be divergent part from the finite part we rewrite this 
integral as follows [^] : 
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iV(p) = [TiM - n^(o) - ^fH%rv(o)] + [rwo) + ^fl^n^(o)]. (25) 

The first derivative term disappears due to the symmetry of the integral. Note that we 
adopt the dimensional regularization with 4 + 2e dimensions to avoid infra-red singularities 
for zero external momentum. With this infra-red regularization the first term has already 
been finite, so that it can be evaluated as the value of the naive continuum limit (a — > 0). 
Thus we obtain, 

iWp) = [n-*'(p) - n?r (o) - ^^n-'-(o)] + [rwo) + ^ p 2 , CT rv(o)] 
= rc n *-(p)] + [n^f ce (o) + ^d p %n^ ce (o)], 

where cant, stands for continuum. In the last equality, we use the fact that the integral with 
zero external momenta is zero in the continuum dimensional regularization. It is noted that 
the "cont. " term is integrated with the dimensional regularization but the second term 
must be integrated on the lattice in 4 + 2e dimensions. In other words, the first term is 
the contribution of the continuum theory and the second term is that of the lattice theory, 
and therefore the latter is named as "lattice". It is also noted that we use 75 which anti- 
commutes with all 7^ in 4 + 2e dimensions for the calculation of 2-point function. Since 
the final result is independent of the infra-red regulator, so that it does not depend on the 
choice of 75. 



a. Evaluation for the Continuum part 

The continuum part leads to the usual transversal form, multiplied by the function 
Fi/R(s,t) which characterizes the domain-wall fermion. The result is 



n^>) = -/£F x 2 (M) 



al n q 



-i(k +Ip) 



(27T 



9 Yf 2 



1 5 , Anfi 2 . 

- + --1E + log — ^ 

e 3 p 2 



2 -T 



(26) 



where e 



n - 4 



and \x is the renormalized point and = T^p) is the transverse function 



of p: T^ u (p) = p 2 5 ilu - p^p u . 
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b. Evaluation for the Lattice part 

The Lattice part U l ®* Uce (p) can be divided into two quantities, a mass term of the gauge 
field and a second derivative term of the IT^, which contains the transverse part. 

We consider the mass term n^ t ce (0) = U M (s,t)5^ u , which has the following form : 

IIm(M) = -£ 2 £ ^tr{[d,S^(q) ■ SMUd^iq) ■ S F (q)] ts } X a" 2 

+ 9%t £ -^^K S f\q) ■ S F (q)}ss x a~\ (27) 

Here no sum over fj, is taken. By the rescaling q — > qa and the fact that the integral is 
infra-red finite in 4 dimensions, we obtain 

U M (s,t) = g 2 a~ 2 j^-0- 4 t, {d,[L,(q)] s , t } (28) 

" g2a ~ 2 L ^ tr {[LM * LM] - [LM ■ LM]} ^ ' (29) 

where 

[LM)]s,t=[d,S- F \q) ■ S F {q)] sU (30) 

and 

[A * B] S)t = A 8>t Bt tS , [A ■ B] Stt = S s , t ]T A S)U B U , S . (31) 

u 

Since the first term is equal to zero because of the Stokes' theorem, the mass term of the 
gauge field finally becomes 

U M (s, t) = -g 2 a- 2 £ ^_tr {[L,(q) * L,{q)\ ~ M?) ' L M\} s ,t ■ 
It is easy to check that this mass term satisfies the following identity: 

L L 

]T n M (M) = E n M (s,t) = o, (32) 

t=-L s=-L 

which comes from the fact that the theory with the axial gauge fixing is still invariant under 
s independent gauge transformations. Thus, for the s independent background gauge field 
A[j,{s,p) = A^(p), no mass term is generated by the fermion 1 loop integral. 
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Since it is difficult to calculate the mass term analytically for general s and t, we evaluate 
it numerically In Fig.|I], the behavior of the mass term n^(s,t) with t = fixed is plotted 
as a function of s at mo = 0.5 and L = 5. The Um{s, 0) has the largest (negative) values at 
s — 0, the place where a chiral zero mode lives. This means that a loop of the chiral zero 
mode mainly contributes to il^(s,t). Furthermore we have checked that Um{s,s) is small 
at s 7^ or ^ L, where only massive modes exist. The behaviour of Um{s,0) is similar to 
the shape of the mass term in 2+1 dimensions f|, which is given in Fig.|^ where D(s, 0) is 
plotted as a function of s at m = 0.5 and L = 5. Here D(s,t) = J2x Fx(s,t) 2 — 2K(s,t) 
and K(s,t) is given in ref. ||. 

Next let us consider the calculation of the second derivative term. Since the transversality 



T) T) 

is hidden behind this term, we can parametrize ^^d 2 p>(T U l ^ tice {0) as follow. 



Pp^ d ^ U laUice {0) = U (a) T ^ (p) + + 11%^ (33) 

The second and the last term in the above equation break the transversality. Especially, 
the last term breaks the Lorenz invariance as in the case of Wilson- Yukawa formulation for 



lattice chiral gauge theories [23|. Although the details of calculation given in Appendix B 
is important, we give only the results of lT a ' 6 ' c ) below, where a and (3 are not equal and no 
sum over them is taken. 

n (a) = -nW + ^£^tr{[(L a (g)-L,(g))^(L /3 (g)-L a (g))] 

-[(Lr(g) " L?»{q)) * (L^(q) ■ L™(q))]}^ 

9 2 ^ p2, Jl_ J_ , 2 x 0.46349 \ 
^)\h/ X \Z V3n + 3tt» ) 

+ ^ x 5/x(^*)^0 + M«V)) ) (34) 

where L c ° n (q) is a value of L^{q) in the naive continuum limit. Suppressing the extra 
dimension indices s, t, becomes 

n (b) = =f £ ^tr {L a (q) ■ L a (q) ■ [L p {q) * 1 - 1 * L p {q)\ ■ L p (q)} 
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+ 4" L Wr ix {La{q) ' [Laiq) * 1 " 1 * L ^ L ^ q) ' + L ^ q) ' 

+ "4" £r (2^ tr ' * 1 - 1 * L °W • L ^ ' ^ 

and IT^ becomes 

n (c) = _ n (b) + 

- |- £ • L a (q) ■ [L a (q) * 1 - l*L«(g)] • 

+ ^tr{[(<9 Q L a (g) • L a (q))*L a (q) 

-2(L a (q) ■ d a L a (q)) -k L a (q) + L a (g) * (d a L a (q) ■ L a (q))]} 
+ [H. C\. (36) 

Here iJ. C. stands for hermitian conjugate. It can be seen easily that IT^ which appears in 
II^'^ satisfies the equation: 

£<J = 0- (37) 
t 

A similar formula can be found for n( c ) after a simple calculation: 

E n 2 = ^ £ (C^^ {^) 3 } = o. (38) 

since the integral of n dimensions without singularities is zero from Stokes' theorem. 
The transverse term satisfies the following equation: 

E n (a) = 2 x ^ (1 + /o,(aV) + 1.147) . (39) 

It is noted that D>) has l/e term. This is an infra-red singularity, which is canceled in the 
total expression of n M „ iS t(p). This point is discussed also in Appendix B. 

c. Total contribution ofH IJjU (p) 

As n™ n *-(p) and Tl l ^ tlce (p) are obtained in the previous paragraph, the total contribution 
of Hfivip) (s,t are suppressed) is now given by 
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n^(p) = [n^-( P )] + [n M + rt a %„( P ) + u^5, u p 2 + uPs^fi] 

= n M + U(p)T^(p) + U^5, uP 2 + U^S^pl, (40) 

where H(p) is the total contribution for the transverse term, which has no 1/e term (see, 
Appendix B). Since it is difficult to calculate eq.fl4"0|) analytically, we numerically evaluate 
each term in eg. (j33|) : a finite part of IT( a ) (A term), lT b ) ( B term ) and IT( C ) ( C term), and 
plot results in Fig.|3] ~ [5| These figures are written as a functions of —L < s < L (L = 5) 
for t = fixed. 

From the above result, we draw the following properties for the structure of the gauge 
field 2 point function at a fermion 1 loop level, as a conclusion of this subsection. 

1. If the gauge field has no dependence of the extra dimension, (i.e. A^s, n) — ► A^{n)), 
the propagator n^(p) becomes transverse, thus gauge invariant, as expected. This 
can be easily seen from a fact that only the transverse term in eq.([E]) is left after 
summation over s,t (using eq. (|32|) , (|37|), (|38|) and (|39"D): 



£iW?) = jSji 2x l TM (41) 



X 



^ + 7£ + (log(^) + 1.147) 



2. IIm, n(p), and n^ c ^ have a peak around s = for fixed t = 0. This shows that 
there is no gauge invariance when the gauge field depends on the extra dimension s, t. 
This fact will be discussed later. 



B. Anomaly in 4+1 dimensions 

The next order of the effective action is a three point function. An important quantity 
which we should evaluate is a divergence of the gauge current, i.e. gauge anomaly. In this 
subsection we shall concentrate on a calculation of the gauge anomaly. An effective action 
for three gauge fields induced by a fermion loop integral is written as 
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xr^ p (q,t,k;s,t,u) x A M (s,g)A,(t,p)A p (w, fc), (42) 
where T fll ,p(q, t, /c; s, t, w) is the three point function of the gauge field. 

(3) 

The anomaly is defined as a variation of S e ^ under the infinitesimal gauge transformation 
5A p (k,u) = ^-q p 6(k,u): 



siu 7r/a (27TJ v y V J 

xik p Tp Up (q, t, k; s, t, u) x A p (g, s)A,(p, t)9(k, u), (43) 
and a relation between SS^j and the current divergence is given by, 

SSi% = V -^<L4 P (A;, M )- 



^ r /a J^L 6A (k u) 5S $f 

rn/a d n k 



More explicitly, 

/7r/a (7™ <7 n r> , , 

-^--^(27r) n (5 (n) (g + p + fc) 
-Tr/a (2tt)™ (2tt)™ v ; 

xk p T pup (q,p,k;s,t,u) x A p (q, s)A u (p,t). (44) 
Since r pl , p has the logarithmic divergence, we evaluate it as follows. 
kpTfiup\Q,P, k] s, t, u) 

= k p [T pi/p (q, t, k; s, t, u) - 1^(0, 0, 0; s, t, u) 

-q a d a T pvp (0, 0, 0; s, t, it) - p^r^O, 0, 0; s, t, it)] 
+A: P [r^ p (0,0,0;s,t, M ) 

+g ff d ff r, MV (0, 0, 0; s, t, u) + p tr 9 ff r AtI/p (0, 0, 0; s, t, «)] . 

We can then evaluate the first term using the continuum form as in the previous subsection. 
Only the parity odd part of T pvp is necessary to get the anomaly, and it becomes 
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kpT ^p^q, p, k; s, t, u) — k p [T pi/ p(q, t, k; s, t, u)] 

+ k p [q a d a Tp U p(0, 0, 0; s, t, u) 

+ Pf7 d a Tp U p(0,0,0;s,t,u)] latUce . (45) 

Here we use a fact that the anomaly has an anti-symmetric structure for Lorentz indices, 

fl,U,p,(T. 

We now perform the calculation of these terms, continuum and lattice ones, in the next 
two paragraphs, and the detail is given in Appendix C. 

d. Evaluation for the Continuum part 

Althouh one can use the dimensional regularization as in the previous subsection to 
calculate the continuum part of the anomaly, we introduce a different regularization here - 

/oo rA 
is replaced by / in n = 4 (A = — ). Accordingly, using = F x (s,t)F x (t,u)F x (u, s), 
-oo J— a a 

we can write 

k p [V p Vp (q,p,k] s,t,u)] cont - 

1 M d 4 l f „ -il n +P 1 ) , „ -iil -h )] 

+(((-> f-p) + ((->(- 5)] 

Here we use the Ward-Takahashi identity as usual in the last step. It is noted that we can 
not make this integral zero using the shift I —> I + q etc, since the integral with respect to I 
has boundaries, —A or A. 

Using the following integration formula: 

rl rA-xp dH 8 1 1 

Jo dx L- xp J^W a pTJ^ = ^ {A ^ oo) > (46) 

where h 2 (x) is a positive function ofx, we finally obtain 



t 

[k P Tp Up (q,p, k; s, t, u)] = ^2 &xF x e ixv P a%Va- 

(47) 
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e. Evaluation for the Lattice part 

The lattice part of the anomaly is written as 



k p \q a d a Yp Vp + p a d a T p Up ] lattlce (0, 0, 0; s, t, u) 

= ^^u P aq P Pa\p(s, t, it), (48) 



where Ag iP (s, t, u) = A q (s, t, u) — A p (s, t, u) and 



1 d 

A q (s, t, u) = —e^pe— r^(0, 0, 0; s, i, u). 



An explicit form of A qtP (s,t,u) is given as 

1 If* dH 

A Q:P (s, t, u) = -g 3 e^ pa - J ^ j—y^tr {(L„(l) ■ L a (l)) * L v (l) * L p (l)} 

-d CT tr{L M (0*M0*M0) (49) 

where the extra dimension indices s,t,u are implicit in the above equation and [Lp] st (l) = 
[dpSp 1 ^) ■ Sp(l)]st is used again. 

As seen in the calculation of the two point function, there are two useful formulae: 

A q , p {s, t, u) = -24^2 Z 5 x F x( s ' fyiwtxrQpP* ( 5 °) 

and 

J2A q , p (s,t,u) = ^E^xK^W^P- ( 51 ) 

The detailed proof is given also in Appendix C. 
/. Results 

A total contribution of the parity-odd term, the gauge anomaly, becomes 

odd 



kpT fxvpijli Pi kj Si t-, U) ^pupaQpPa 



X 



ig 



3 



24tt» ^ 



(52) 



We summarize important properties of the anomaly as follows. 
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1. Note that a summation over u makes this anomaly zero due to eg. (|50|) . This comes from 
the fact that the model with the axial gauge fixing is invariant under s independent 
gauge transformations. 

2. Because of eq . (|5l|) , a summation over s,t makes this anomaly equal to 

k p V nup{Q? Vi k'i S, t, u) \gnomaly = 3 2 ^HVpcrQ.pPcr ^X^X U ) ' (^3) 
st O 71 " X 

Physically this is the gauge anomaly at u for s independent background gauge field 

1 

A^s.n) = A^in). It is noted that the coefficient of eq. (|53|) , — -, is equal to the 
coefficient of the covariant anomaly in 4 dimensions, not the one of the consistent 
anomaly in 4 dimensions. This does not contradict with the fact that the above 
anomaly is derived from the variation of the effective action, since the effective action 
is defined in 4 + 1 dimensions, not in 4 dimensions. 

3. Without summations it is difficult to calculate the anomaly analytically. In- 
stead we evaluate it numerically, and results are given in Figs.[§-|8|. We plot 
[ — 2 5 x F^(s, t, u) + A q:P (s,t,u)] as a function of u at m = 0.5 and L = 10, 
for fixed s = t = in Fig.^ for s = and t = 2 in Fig.[7|, and for s = and t = 8 
in Fig|| These figures tell us the following. If two gauge fields A M (s) and A v (t) are 
on the same 4 dimensional subspace (s = t = 0), the anomaly has large contribution 
around u = s = t. If t differs a little from s ( s = and t — 2 ), the anomaly has 
non-zero contribution around u = s or u = t but peak heights become less. If t is far 
away from s (s = and t = 8 ), the anomaly almost vanishes at all u. For comparison 
we calculate the anomaly in 2+1 dimensions, which is given by [f| 

where 

C(s, t) = Y, SxF x (s, tf - 2T CS ( S) t). (55) 
x 
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In Fig.^ C(s,t) is plotted as a function of s with t = fixed at mo=0.5 and L = 5. 
Behaviors of the gauge anomalies are similar both in 4+1 dimensions and in 2+1 
dimensions: they have large contribution around s = or u = 0. 

Finally it should be mentioned that a parity-even part of the variation of the 3 point 
function under the gauge transformation vanishes locally for s independent background 
gauge fields. 



C. The 4-point function 



A calculation of a 4-point function is simpler than the previous ones since there is no 
derivative term. Using the same logic as before we obtain 



r 2L(^ P, k, r; s, t, u, v)= frjfj (q, p, k, r; s, t, u, v) 



cont. 



+ 



r£> {0,0,0,0; s,t,u,v) 



lattice 



Results are very similar to the 2- or 3-point functions; After the summation over v 
the gauge variation of the 4-point function is zero and after the summation over s, t, u it 
reproduces the continuum value. It shows that the gauge invariance is correct only in the 
case that the gauge field has no dependence on the extra dimension. 



IV. PHYSICAL IMPLICATION 

We consider the 4+1 dimensional theory so far. Properties of the effective action in 
4+1 dimensions are more or less similar to those in 2+1 dimensions except that the 2-point 
function of the gauge fields contains a mass term and a Lorentz non-covariant term as well 
as a transverse term and a longitudinal term, and that the gauge anomaly is proportional 
to the charge cubed. As in the 2 dimensional case the fermion number violation can be 
incorporated as follows. First let us consider a fermion number current, whose expectation 
value for background gauge fields is defined by 

W'-*» = 5£fe (56) 
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where the index g in the current explicitly shows the charge of the fermion, and the gauge 
current is equal to gJ2(s, x). A divergence of the fermion number current, the fermion num- 
ber anomaly, is proportional to g 2 while the gauge anomaly is proportional to g 3 . Therefore 
even for anomaly-free model such as the standard model, where a sum of charge cubed 
vanishes, the fermion number violation, which is proportional to a sum of charge squared, 
can remain non-zero in the domain-wall fermion formulation. This is the most important 
consequence of our results: The domain-wall fermion formulation may allow us to simulate 
the non-perturbative dynamics of the baryon number violation of the standard electro-weak 
theory. Although our results are obtained for a U(l) gauge theory, it is not so difficult to 
obtain similar results for non-abelian gauge models. 

The gauge non-invariant terms remain non-zero in 4(+l) dimensions even for anomaly- 
free cases as in 2(+l) dimensions. In the next section we will briefly mention a way to avoid 
these terms. The s-independent gauge fields again lead to the vector-like theory as in the 2 
dimensional case |J. 

V. CONCLUSION 

In this paper we have investigated the domain- wall chiral fermion formulation in 4+1 
dimensions with a lattice perturbation theory. We have calculated in detail the 2-point 
function and the gauge anomaly in 4 dimensional U(l) chiral gauge theory formulated in 
4+1 dimensions. 

The most important conclusion drawn from our perturbative analysis is that the baryon 
number violation of the standard model may be incorporated by the domain-wall fermion 
formulation. 

However, as stressed in the previous section, the domain-wall fermion formulation can 
not maintain the transversality, and hence gauge invariance, even for anomaly-free cases. 

A solution to the gauge invariance has been proposed by Narayanan and Neuberger: 
They take the size of the extra dimension strictly infinite, instead of periodic box, and make 
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gauge fields independent of the extra dimension. (They also take the "temporal gauge" 
from the beginning.) Since no anti-domain wall caused previously by the periodicity of the 
extra dimension exists, the theory is expected to remain chiral even for such gauge fields. 
Indeed, using our result of the effective action in 2+1 dimensions it has been shown that the 
Narayanan-Neuberger formulation, called "overlap formula" , gives gauge invariant effective 



theory except for the gauge anomaly in 2+1 dimensions [ 14| . The U(l) chiral gauge theory 
via their formulation in 4+1 dimensions is now under investigation since our calculation in 
this paper can be easily extended to it. As mentioned in the introduction, however, it is 
pointed out recently that the overlap formula can not reproduce the desired chiral zero mode 
in the presence of the rough gauge dynamics ( gauge degree of freedom ) which appears due 
to the violation of the gauge invariance at infinities. This problem is known to exist also for 
the wave-guide model. 

Except the problem of the gauge invariance above, the domain-wall fermion formulation 
in 4+1 dimensions works well for smooth back-ground gauge fields. A remaining question 
of the domain-wall fermion formulation is whether the chiral zero mode can survive in the 
presence of the dynamical gauge fields. A numerical investigation has been performed in 2+1 
dimensions but has failed to give a definite conclusion on an existence of the zero mode in 
the symmetric phase [13|. Further numerical study for a 4+1 dimensional model is needed. 
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APPENDIX A 

In this appendix, we calculate the function Gl/r(p) ■ The fundamental equations are 

p 2 G L (p) s , t + {M^MG L (p)) s , t = 5 Btt , 

p 2 G R (p) s , t + (MM^G R (p)) s , t = 5 Btt , (57) 

where X)«=i sin 2 ^) is written as p 2 . We examine the upper equation in detail. Explict form 
of this equation is 

(p 2 + A(s))G, tt (p) - B{s + l)G(p). +1 , t - B(s)G(p)._ lft = S a , t , (58) 

where 

A(s) = l + B(s) 2 , (59) 
B(s) = a(s) - r -F(p) 

= l-m( S )-^(cos(^)-l). (60) 

Since m(s) in B(s) changes the value at s — 0, L, eq.(^) should be separated into two cases, 
one for < s < L — 1 and the other for — L — 1 < s < —1. For convenience, G(p) + is defined 
as G(p) in the range of0<s<L — 1 and G(p)~ in — L — 1 < s < — 1. 
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We first focus our attention on G + (p). In the range < s < L — 1, using B(s) = 
B(s + 1) = 1 — m — 7}F(p) = a + , eq.(|58D is rewritten as follows: 

(p 2 + 1 + a 2 + )G+ (p) - a+(G(p)+. 1>t + G(p)+ 1|t ) = (61) 

The solution of this equation is expressed as a sum of homogeneous general solutions and 
an inhomogeneous special one. Homogeneous general solutions with two unknown functions 
g+ ■* (t) and g+ (t) are 



9+ 



1 1 + p 2 

where cosh(a + (p)) = ~(a + + 



2 V a+ 7 

Next the inhomogeneous solution is calculated with the Fourier transformation as (l n 



2nn/L) 



I L - 1 e lln\s-t\ _J L-l g\g- t \+l r?\ s - t \ + l 



2La+ ^ cosh(a+(j»)) - cos(Z n ) 2a + Lsinh(a + (p)) ^ L Z - e Q +W Z - e- Q +(p)' 

cosh(a+(p)(|s-*| -L/2)) 
2a + L sinh(a + (p)) sinh(a + (p)L) ' 

where the following formula is used in the last step: 



(63) 



z£— ^ (W ">- (64) 



Therefore, the solution G(p)f t is given as 



+ cosh(a + (p)(|s-t| — L)) 
2a + Lsinh(a + (p)) sinh(a + (p)L) 



In same way, the solution G(p) st is obtained as 



G(p)- t = g w (t)e a ^ s + g { '\t)e- a -^ s 

+ cosh(q_(p)(|s-t| — L)) 
2a_L sinh(o;_(p)) sinh(a;_(p)L) ' 



1 1 + p 2 

where a_ = 1 + m — cosh(o;_(p)) = -( a - H )• 
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The unknown functions g±\t) and g±\t) are determined by the four boundary condi- 
tions, which are obtained by considering eq.(^) at s = and L: 

G(p)+_ 1>( = G(p)~ = _ lt , 
a+G(p)f=o,t = a-G(p)j =0jt , 
G{p)t=L-x,t = G (p)7=-L-i,t, 

a+G(p)t=L,t = a ~G{p)~ = _ L)t . (67) 



Explicitly, 



e 



3 o+ 



-a+ 



e 

a + a + 

a+e -a+L a+e a+L 
e -OL+(L-\) e o+(L-l) 



— e 



— e 



-a_ 



-a_L 



-a_ 



-a_e — a_e 

. e -o-(H-l) _ e a_(i+l) 



sl 2) (t) 

9^ it) 



b(t) 



\ 



X_ cosh(«_(|l + t| - L)) 

-X+cosh(a+(|l + t| — L)) 
a_X_ cosh(a_(|t| — L)) 

—a + X + cosh(a + (|t| — L)) 
a_X_ cosh(a_(| — L — t\ — L)) 

—a + X + cosh(a + (|L — t\ — L)) 
X- cosh(a_(| — L—l — t\ — L)) 

-X+ cosh(a+(|L — 1 — t\ — L)) 
where X^ 1 = 2a± sinh(a±) sinh(a±L). Although solution to these equations is very compli- 



V 



(68) 



cated in general, it becomes simpler in the limit of L — > oo. For t ~ 0, 



e 



-a+ 



\ 



a + a + — a_ 
a + e a + L 
e a +( L -^ 



— a_ 

-a_e a ~ L 

_ e a_(L+l) 
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cosh(a_(|l + t| — L)) — X + cosh(a + (|l + t\ — L)) ^ 
a_X_ cosh(a_(|t| — L)) — a + X + cosh(a + (|t| — L)) 






Since g± — > in the limit of L — > oo, g± is easily obtained from this matrix. It is noted 
that at t ~ L, — > in the limit of L — > oo. This shows that Gi(p) has contributions 
from two chiral zero modes which live at s = or L. 



APPENDIX B 

In this Appendix B, we show in detail the calculation of the 2-point function 11^ (p) in 
4+1 dimensions. 

Let us set \i = a and v = (3 (a ^ (3) in the parametrization of eq.(^3|). The extra 
dimension indices s,t are suppressed below. For this choice of \x and v eq.(|33]) becomes 



5^^ T BS«»(o) = nW(-p a p^ 



It means that 



(69) 



(No sum over a and (3.) When we set \i = v = a, eq.(|33l) becomes 

^a;, ff nr(o) = (nw + n( & ))(p 2 - P D + (nw + 

Therefore relations between IT 6 )'^ and n MI/ (p) are given by 

n( 6 ) = i^n aa + ^XT ce (o), 
n< c > = -^n^* iice (o) -n (fe) . 

2 



(70) 



(71) 
(72) 



The most simple term in eq.([33|) is IT 6 ), which explicitly given by 
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8 «/ — 7t (2vr)" 



tr {[<9^L Q (g) *L Q (g) 



- 2 [d^L^g) * ^L a (g)] + \L a {q) -k <9jL a (g) 
+ 2 ^L aQ (g) * Lg(g)] - 2 [L Qa (g) * L^g)] 

- 2 [fyL a (g) * a 2^(g)] + 2 [L a (g) * ^L^g)] } 



(73) 



where [L M ] st (g) = [c^S^g) • S F (q)] st and [L w ] rt (g) = [d fl S F 1 (q) ■ dpS F {q)]irt- We obtain 
eq.(|35D by substituting the following relation into the last equation. 



L m (q) = d'S F \q) ■ S F (q) + L,(q) ■ L,(q). 



(74) 



Here we use the fact that the following n- dimensional integral is zero by the Stokes' theorem: 

£^tr^[L o (g)*L o (g)] = 0. (75) 

Il( c ) is also obtained in the same way. Note that both 14^ and IT C ) have no infra-red 
singularity. 

Next we consider the most interesting term U^ a \ For convenience eq.(|69D is rewritten 
with eq.([n]) as 



n w = -n( & ) + ^n aQ (o). 



(76) 



By this identity we can concentrate on the last term for the calculation of U^ a \ The last 
term becomes 

— g 2 r n d n q 



-7T 2tt 



g 2 f w d n q 



-tr { [d 2 L a (q) * L a (q)] - 2 [d p L a (q) * d p L a (q)] + [L a (q) * d 2 L a (q)] } 



-7T (2 



7T 



-trdp {[dpL a (q) -k L a (q)} - [L a (q) -k dpL a (q)]} 
-4 x tr {[dpL a (q) -k d p L a (q)]} . 



(77) 



To derive the last equation we use the formula (|73]) again. By the Stokes' theorem in n 
dimensions, this integral simply becomes 

g 2 n d n q 



n(°) + n (fc) 



2 J — 7r (2tt)' 



-tr{[^L a (g) *dpL a (q)]} 



(78) 
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Since this integral has a logarithmic infrared singularity we evaluate it in the following way: 

_|_ ]j( b ) = _ ttO)| co ™*- _|_ tt<»| co ™*- _)_ ]j( b ) ) (79) 

where n( a ) |cont. means ^he integrand in lT a ) is replaced by the continuum one: 

-g 2 1^ 2/ \ r d n q \ . — ih . — ^ . — yt[ . —ify 1 



n( a) |CO n, = ^ M) V f2(s t) r _o_ tr • zw -w z# . 

1 8 1 * } 2^ x^'VJ-* (2tt)" \ 9 2 ? 2 ? 2 ? 2 



8-2 ^ xV ' 7 7-. (2ir) n [ (g 2 ) 4 (g 2 ) 2 J ' 
and ]j( b )\ cont - = o. Using a property of the Gauss's error function, we obtain 

n<.>r<- = VE**) fef'W^ i >.)" + (5 - 

where < 1 > y is defined by jl^g dxle~ x2 . The first term in the last step has a 1/e divergence. 
Therefore we evaluate this integral as follows: 



80) 



y) 



poo rl 

/ dyy l - n '\< 1 > y ) n = / dyy-\< 1 > 
Jo Jo 

+ J™ dyy- 1 {(< 1 > y f - v^F 4 } + v^F 4 dyy l - n 



-7T 2 

— + 0.46349. 



Here we set n = 4 for the term without infra-red divergence. 



Finally we reach the result eq. (|34]) . By considering that the factor a n jl n must be mul- 
tiplied in front of the n dimensional integral, it can be easily understood that log(a 2 jl 2 ) 
appears in eq.(^) |2"2fl . This is one of the most important points among what the final 
result tells us: The e -1 of the lattice contribution cancels that of the continuum one and Jx 
dependence also disappears, so that only the ultra-violet divergence logp 2 a 2 remains in the 
final result. We can also see a similar kind of the cancellation of infra-red divergences in 
Ref. M\. 



APPENDIX C 



We investigate the anomaly in this Appendix. Let us start with the following definition 
of the three point function: 

28 



1 



dH 



3 J-. (2tt) 4 



rtr 



{[dpS-^l-qfflSil)]^ 
d u S~\l + p/2)S(l + p)\ *\d p S~ 1 (l + (p-q)/2)S(l-q)} } 

J tu L J MS J 

+ (/^/-j9) + (/^/-g). 



Substituting the above into eq. fl48|) , we obtain eq.(|49|), which is given once more below: 



dH 



3 J-tt (2tt) 4 



(81) 



x3tr {[L^l) ■ L a (l)] st * L u (l) tu * L p (l) us } - <9 CT tr {L p (l) st * L v (l) tu * L p (l) us } . 

The dependence on the extra dimension s,t,u of this anomaly is calculated numerically. 
On other hand, we can calculate it analytically in special cases - sum over s,t and sum over 
u. For considering these cases, let us prove two identities, eq.(|50|) and (|51|). First we show 
that 




Aq. p (s, t, u) 



^ 9 e pvpa 



- dH 1 



-7T (2vr) 4 3 

-Ad a tT{[L,(l).L u (l).L p (l)] uu } 
+2d a tr{[L p (l)-L u (l)] st *L p (l) ts } 

using the following line of identities to the first term of eq.(|8TD (s,t,u are suppressed): 



(82) 



s,t 



£« 



3ti{[{L^l)-L a {l))*L v {l)*L p {l)}} 



3tr{[L p (0-^(0- L a {l) ■ L V {1)}} 
3tr {[L p (l)-L a (l)}* [L u (l) ■ L p (l)}} 

-3tr{[^- x (0 • d p S(l) ■ drS- 1 ® ■ d u S(l)}} 
+3tr{[c^- 1 (Z) • d„S{l)\ * [drS- 1 ® ■ d p S{l)}} 

Since integrands in the right hand side of eq. (^2|) become total derivatives, the integrals 
can be performed with the Stokes' theorem. Using the equation 
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[LMst = *7„ x -4- E F x& t)P x , (* - 0) 
and introducing the infra-red cut-off e, the total derivative is evaluated as 

^^{^(0*^(0*^(0} 





d 4 Z 




U (2tt) 4 


1 3 


/ >7r d 3 Z 


— 4!^ e A t ^ CT 


l-n (27r) 4 



2 x ii4! e V ^ rf Eoja 
In the last equality we used the following formula: 



Z CT =7r/2 

— I 



X 167T 2 
Z CT =e x 



Eq.(||) and (0) lead to identities eq.® and (0). 

We shall turn to the proof of e pupa d a Y ^piQ, 0, 0; s, t, u)\ cont - = in eq. , which is 
nontrivial in our regularization, using the same technique above. The calculation is almost 
the same as that for eq. flST|) except that the integral — > J_ A and L^l) is that of the 
continuum: 

/ivp(0, 0, 0; s, t, u)\ con ' 

= W~ E *J tSi 3 ^ • M0 • M0 • M0) - d a ti {L,{1) ■ L V {1) ■ L p (l)} 

Therefore we should calculate the following integral instead of eq. ( p3|) . 

A r/V /„ ' S 



l„=e 



-a (2tt) 4 (Z 2 ) 2 
La (2tt) 4 7o 



r A/e r 1 dxxe ^' 

J-A/e (2vr) 4 Jo 



1 1 , 

0. 



167T 2 ' 167T 2 ' 

This means that the derivative term in the continuum, e Ati/pcr 9 (r r Mi/ p(0, 0, 0; s, t, u)\ cont , is zero 
in this regularization. 
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FIGURES 



Mass term 

Source at t=0 , amO = 0.5 
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s 

FIG. 1. The coefficient of the 2 point function Ilfcf (s,i) as a function of s with t = fixed at 
mo = 0.5 and L = 20. 
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D(s,t) 

Source at t=0 , amO = 0.5 



2.00 



1.00 
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-2.00 ' 1 1 1 1 1 1 1 1 1 1 
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s 

FIG. 2. The coefficient of the even term D(s, t) in 2+1 dimensions as a function of s with t = 
fixed at mo = 0.5 and L = 5. 
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A term 

Source at t=0 , amO = 0.5 
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-5.0 -4.0 -3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 4.0 

s 

(a) 



5.0 



FIG. 3. The coefficient of the 2 point function ir st ' as a function of s with t = fixed at 



mo = 0.5 and L = 5. Here divergent contributions are removed from n^. 
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B term 
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FIG. 4. The coefficient of the 2 point function 11^ as a function of s with t = fixed at 
niQ = 0.5 and L = 5. 
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FIG. 5. The coefficient of the 2 point function as a function of s with t = fixed at 
niQ = 0.5 and L = 5. 
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Anomaly 

Source at s=0,t=0, am0=0.5 
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FIG. 6. The coefficient of the anomaly C(s, t, u) in 4+1 dimensions as a function of u for s, t = 
fixed at mo = 0.5 and L = 10. 
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Anomaly 

Source at s=0,t=2, am0=0.5 
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FIG. 7. The coefficient of the anomaly C(s, t, u) in 4+1 dimensions as a function of u for s = 
and £ = 2 fixed at mo = 0.5 and L = 10. 
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Anomaly 

Source at s=0,t=8, am0=0.5 
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FIG. 8. The coefficient of the anomaly C(s, t, u) in 4+1 dimensions as a function of u for s = 
and t = 8 fixed at mo = 0.5 and L = 10. 
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FIG. 9. The coefficient of the anomaly C(s,t) in 2+1 dimensions as a function of s for t = 
fixed at mo = 0.5 and L = 5. 
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